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Abstract 

Measurements performed at the Tevatron of both the like-sign dimuon charge asymmetry 
in Bd jS -meson samples and the mixing-induced CP asymmetry in B s — > J/ip(ft depart from 
their standard model (SM) predictions. This could be an indication for new CP phases 
in AB = 2 transitions, preferentially in B s -B s mixing. The experimental situation, 
however, remained inconclusive, as it favored values of the element T\ 2 of the decay 
matrix in the i? s -meson system that are notably different from the SM expectation, 
suggesting the presence of new physics in the AB = 1 sector as well. The very recent 
LHCb measurements of B s — > J/ip(f> and B s — > J/tpfo, which do not find any evidence 
for a new-physics phase in the element M[ 2 of the mass matrix, point into this direction 
as well. In this article, we explore the potential size of non-standard effects in 
stemming from dimension-six operators with flavor content (s6)(fr). We show that since 
the existing constraints imposed by tree- and loop-level mediated .B^-meson decays are 
quite loose, absorptive new physics of this type can in part explain the current data. 
Our model-independent conclusions are finally contrasted with explicit analyses of the 
new-physics effects in B s —B s mixing that can arise from leptoquarks or Z' bosons. 
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The phenomenon of neutral 5 s -meson mixing is encoded in the off-diagonal elements Mf 2 and 
r^ 2 of the mass and decay rate matrix. These two complex parameters can be fully determined 
by measuring the mass difference AM, = M^ — M[, the CP-violating phase </>j/^ = — 2/3 s , the 
decay width difference Ar s = T S L — T S H , and the CP asymmetry a s j s in flavor-specific decays. 
The combined Tevatron and LHC determination of the mass difference reads [TJ [2] 

AM S = (17.73 ± 0.05) ps _1 , (1) 

and agrees well with the corresponding standard model (SM) prediction [3] 

(AM S ) SM = (17.3 ± 2.6) ps" 1 . (2) 

Here the quoted errors correspond to 68% confidence level (CL) ranges. 

The phase difference (fij/^ between the B s mixing and the b — > sec decay amplitude and 
the width difference AT S can be simultaneously determined from an analysis of the flavor- 
tagged time-dependent decay B s — > J/tfj<j). The latest HFAG averagd^ combines the results of 

1 Notice that only results published or accepted in a refereed journal by March 15, 2011 have been included 
in the average. 
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CDF [HE] and D0 [6], which are based on 1.35 fb" 1 and 2.8 fb" 1 of pp data at y/s = 1.96 TeV, 
respectively. At 68% CL they find [7J 

= (" 44 -2D° , = (0.154ttg5S) ps" 1 , (3) 

as well as a symmetric solution located at 180° — and — Ar s . The corresponding numbers 
in the SM are Efl 



(0j/<W>)sm = arg 



{v c %v cb f 



-2.1 ±0.1)°, (Ar s ) SM = (0.087 ±0.021) ps" 1 . (4) 



They deviate from (|3| by 2.3a and 0.9a, respectively. Both the CDF pi)] and D0 [TT] col- 
laboration have upgraded their B s — > J/ip(J) analysis to integrated luminosities of 5.2 fb -1 and 
8.0 fb" 1 . Focusing on the solution with positive Ar s , the weighted average of these measure- 
ments reads <f> a Jf ^ = (-32 ± 18)° and Ar, = (0.097 ± 0.032) ps" 1 at 68% CL. Compared to 
(j3]) these results show a better agreement with the SM expectations Q, because the tensions 
in 4> s j/^ and Ar s are reduced to 1.7a and 0.3a. These measurements were very recently 
superseded by a LHCb analysis based on 337 pb" 1 of data collected at y/s = 7 TeV [T2] . Com- 
bining the B s — > J/ip<p [33] and B s — > J/ipfo [13] channels the LHCb collaboration finds no 
evidence for a new-physics phase in B s -B s mixing and a value for Ar s that is 1.0a above the 
SM expectation Q. The corresponding central values and 68% CL error ranges are 

S W = (1.7 ±10.0)° , Ar s = (0.123 ±0.030) ps" 1 . (5) 

The CP asymmetry in flavor-specific decays a s ^ s can be extracted from a measurement of 
the like- sign dimuon charge asymmetry A\ h) which involves a sample that is almost evenly 
composed of B& and B s mesons. Performing a weighted average of the results by CDF [15] 
and D0 [T6l IT?] , which are based on 1.6 fb -1 and 6.1 fb" 1 of integrated luminosity, leads to 

4 L = (-8.5±2.8)-10- 3 , (6) 

at 68% CL. Utilizing the SM predictions for the individual flavor-specific CP asymmetries [3] 

(4Jsm = -(4.1 ± 0.6) • 10" 4 , (a%) SM - ( I .9 ± 0.3 i • 10~ r> . (7) 

one obtains 



(4l)sm = (0.506 ± 0.043) (aJ s ) SM + (0.494 ± 0.043) (a* J SM = (-2.0 ± 0.4) • 10" 4 , (8) 

which differs from ^ by 3.0a. Using the measured value of the CP asymmetry in flavor- 
specific B d decays [7J, 

aJ s = (-4.7±4.6)-10- 3 , (9) 
and (|6j) in (|8j) , one can also directly derive a value of a s ^ s . One arrives at 

a s fs = (-1.2 ±0.7) • 10" 2 , (10) 
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A recent concise review of the theoretical uncertainties plaguing (Ar s )sM has been given in [5], 
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which compared to ^ represents a discrepancy of 1.7a. Recently, the D0 collaboration has 
updated its analysis of the like-sign dimuon charge asymmetry. Employing 9.0 fb -1 of data, 
they now find [18] 

4 L = (-7.87±1.72 stat .±0.93 syst .)-10- 3 . (11) 
This new number translates into 

(4l)sm = (0.594 ± 0.022) (a d fs ) SM + (0.406 ± 0.022) (a s fs ) SM = (-2.4 ± 0.4) • 10~ 4 , (12) 

and 

a s fs = (-1.3 ±0.8) ■ 10~ 2 , (13) 

corresponding to a tension with a statistical significance of 3.9<r and 1.5cr, respectively. 

In view of the observed departures from the SM predictions, it is natural to ask what kind 
of new physics is able to simultaneously explain the measured values of AM S , 0}/^, AI\,, 
and a s j- s . Assuming that the found anomalies are solely due to non-standard effects in the B s - 
meson system]^] this question can be addressed in a model- independent way by parametrizing 
the off-diagonal elements of the mass and decay rate matrix as follows 

M S X2 = (M S 12 ) SM + (M* 2 ) NP = (M^sm Rm e l<t>M , 

(14) 

Ti2 = (ri 2 )sM + (r^NP = (T S 12 ) SM Rr e 1 ^ . 

In the presence of generic new physics, the 5 s -meson observables of interest are then given to 
leading power in |r| 2 |/|M* 2 | by 



AM S = (AM S ) S M Rm , <Pj/^ = {<Pj/^)sm + <Pm , 
Ar s = (Ar s ) SM Rr COS ^sm + 4>m - 0r) w (Ars)gM ^ ^ ( _ } 

COS 0| M 

sin (0| M + <j) M -(h) / 5 \ sin (<f> M - r ) 

O-fs ~ { a fs)SM TT- —— s (O/JSM TT" ~ s • 

Um sm( PsM n M (PsM 

Notice that in the case of Ar s and a s ^ s the final results have been obtained by an expansion 
in </>g M = arg (— (Mf 2 ) SM /(r^ 2 ) SM ). Although this is an excellent approximation, given that 
0g M = (0.22 ± 0.06)° [S], we will employ the exact analytic expressions in our numerical 
analysis. 



The four real parameters Rm,t and 0M,r entering (14) can be constrained by confronting 
the measured values of AM S , <\> s ji^ Ar s , and aj s with their SM predictions. In the following, 



3 Hints for deviations also exist in the i^-meson system, in the purely leptonic decays B — > rv T , and in 
CP violation in the neutral kaon system [5J [THl HOI HI] • We believe that in these cases the disagreements are 
not large enough to exclude the possibility of statistical and/or systematic origins of the discrepancies. The 
recent next-to-next-to-leading (NNLO) order calculations of €k [HI [23] , which find theoretical errors that are 
larger than the ones reported at next-to-leading order (NLO), support this suspicion for what concerns the 
neutral kaon sector. In our work, we hence focus on new-physics effects in B s -B s oscillations, ignoring possible 
non-standard effects in other sectors. The main points we want to make do not depend on this restriction. 
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Figure 1: Left (Right): Constraints on the parameter R M and 4>m (Rr an d </>r) from 
the global fit to the data set Dl in scenario SI (S2). For the individual constraints 
the colored areas represent 68% CL regions (dofs =1), while for the combined fit the 
light red area shows the 95% probability region (dofs = 2). In both panels the SM 
prediction (best-fit solution) is indicated by a dot (cross). 



we will analyse two different data setsQ The first set (Dl) consists out of ([T]), ([3]), and (10), 
while in the second case (D2), we rely on <Q, Q, and We be gin our analysis by asking 
how well the SM hypothesis describes the two sets of data. In the case of Dl, our global fit 
returns x 2 = 9-1 corresponding to 1.9cr (2.6a) for 4 (2) degrees of freedom (dofs), while for 
D2 we obtain x 2 — 3.4, which translates into 0.7cr (1.3er) if 4 (2) dofs are considered. One 
observes two features. First, the quality of the fit improves notable when going from Dl to 
D2, which is due to the fact that the value of (frj/^ in the new data set is fully consistent with 
the SM expectation and, second, even the new data is not in perfect agreement with the SM 
hypothesis as a result of the observed large negative value of a s j s (or equivalently A^ L ) and, to 
a lesser extend, the latest measurement of Ar s , which resides above the SM prediction. 

It is also instructive to determine the best-fit solution of the full four-parameter fit to the 
data. In the case of Dl, we find that (Rm, (Pm, Rr, 4>r) — (1.02, —44°, 3.1, 13°) leads to a vanish- 
ing x 2 , while for D2 the x 2 function is minimized for (Rm, <Pm, Rr, 0r) = (1-01, 3.8°, 3.0, 65°). 
After marginalization the corresponding symmetrized 68% CL parameter ranges are 

R M = 1.05 ± 0.16 , d> M = (-46 ± 19)° , 

(16) 

i?r = 3.3±1.5, r = (7±3O)°, 

For simplicity we ignore the mirror solutions at 180° — (f>ji^ an d — Ar s in our discussion hereafter. 
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Figure 2: Left (Right): Constraints on the parameter R M and <p M (R r and 0r) from 
the global fit to the data set D2 in scenario SI (S2). For the individual constraints the 
colored areas represent 68% CL regions (dofs = 1), while for the combined fit the red 
(light red) area shows the 68% (95%) probability region (dofs = 2). In both panels the 
SM prediction (best-fit solution) is indicated by a dot (cross). 



and 

R M = 1.05 ±0.16, <j) M = (1.5 ± 10.0)° , 

(17) 

i?r = 3.4±1.7, r = (58±23)°. 

Focusing on R r and 4>g, we see that the fit to both data sets prefers large corrections to r* 2 , 
either without (Dl) or with (D2) a large weak phase. Based on older data sets this observation 
has been made before in [Ml E5j [26], E7] and triggered a considerable amount of theoretical 
investigations PBIEDIESIEIIISIESESES]. The first publication that, to the best of our 
knowledge, argued convincingly for the possibility to have sizable new-physics effects in rf 2 is 
the work [36] . 

To further elucidate the latter point, we analyse two orthogonal hypothesis of new physics 
in B s -B s oscillations, namely a scenario with (M* 2 )np 7^ and (r^NP = and a scenario 
with (M[ 2 )np = and (rf 2 )NP 7^ 0. We refer to this two different scenarios as SI and S2, 
respectively. The left (right) panel in Figure [I] shows the results of our global fit for scenario 
SI (S2) in the Rm cos 4>m-Rm sin 0m (-Rr cos0r~-Rr sin0 r ) plane based on the data set Dl. 
By inspection of the left plot, we infer that in the scenario SI the preferred 95% CL parameter 
region, is essentially determined by the interplay of AM S and (j> s ji^- The best-fit point is 
located at {Rm, 0m) — (1.01,-41°). It has x 2 /d°fs = 3.1/2 corresponding to 1.3crJ^] which is 
in almost equal parts due to Ar s and aj s . This moderate tension is also clearly visible in the 

5 This result agrees well with the global analysis performed in [5], but disagrees with the work [33], which 
claims that the scenario SI is disfavored by more than 3cr. 
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figure. Turning our attention to the right panel, we observe that in the case of the hypothesis 
S2 only a small region of parameter space has a CL that is better than 95%. The best-fit 
solution resides at (Rr, <pr) = (3.0, 56°) and has x 2 /dofs = 5.5/2 corresponding to 1.9a. These 
numbers imply that the data set Dl prefers the hypothesis SI over S2. 

By repeating the above statistical analysis for the new data set D2, we obtain the results 
displayed in Figure [2] From the left panel, we glean that in the scenario SI the regions 
of all individual constraints apart from aj s now overlap. Minimizing the \ 2 function gives 
(Rm,4>m) = (1.03,2.1°) and x 2 /dofs = 3.3/2 corresponding to 1.3cr, which represents only 
a marginal improvement with respect to the SM hypothesis. In consequence, the case for a 
non-zero non-standard contribution to Mf 2 is rather weak. Comparing the right panels of 
Figures [T] and |2j we see that the 95% CL region is not only significantly larger in the scenario 
S2 than before, but that now also a description of the data with a probability of better than 
68% is possible. The best-fit point is located at (Rr,<j>r) = (2.9,61°). In fact, the latter 
parameters lead to an almost perfect fit with x 2 /dofs = 0.2/2 corresponding to 0.1 cr. The new 
data set D2 hence statistically favors the new-physics scenario S2 over the hypothesis SI. 

Given these results, we believe, that one should take the possibility that absorptive new 
physics could be responsible (in part) for the D0 observation of the anomalous dimuon charge 
asymmetry seriously. While in principle any composite operator (sb)f, with / leading to 
an arbitrary flavor neutral final state of at least two fields and total mass below the B s - 
meson mass, can contribute to Y\ 2 i the field content of / is in practice very restricted, since 
B s — >■ / and B d — > X s f decays to most final states involving light particles are severely 
constrained. One notable exception is the subclass of B s - and I^-meson decays to a pair of 
tau leptons [211 E7J [301 EH EU EEJ E7] , to which we will devote our full attention in this article. 
Specifically, we perform a thorough model-independent analysis of the impact of the complete 
set of (sb)(fr) operators on rf 2 , taking into account all existing direct and indirect constraints. 
We find that the loose bounds on most of the Wilson coefficients of the considered operators 
allow for significant enhancements of r{ 2 compared to its SM value, which help to lessen the 
observed tension in the -Brf jS -meson data. Our work clarifies and extends the existing analyses 
[2H [30], ES] • While most of our discussion is based on an effective-field theory and therefore 
completely general, we also consider two explicit models of new physics that can alter rf 2 , 
namely leptoquarks and scenarios with new neutral gauge bosons. 

This article is organized as follows. After introducing important definitions and notations 
in Section [2j we study in Section [3] the direct bounds on the Wilson coefficients of the (s6)(rr) 
operators following from B s — > t + t~ , B — > X s t + t~, and B + — > K + t + t~ . Section[4]is devoted 
to a comprehensive discussion of the indirect constraints arising from h — > s'j, b — > s£ + £~ 
(£ = e,/i), and b — > 577. In Section [5j we calculate the effects of the full set of (sb)(fr) 
operators on r* 2 and analyze their numerical impact. Our model-independent findings are 
contrasted with two explicit SM extensions in Section |6j A summary of our main results 
and our conclusions are presented in Section [7| A series of appendices contains useful details 
concerning technical aspects of our calculations. 
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2 Preliminaries 



In the SM, the effective AB = 1 Lagrangian is given by 

AG 

£ c ff = £qcdxqed (u,d,s,c,b,e,n,T) + — =- V t *V tb ^ C i(^)Qi , ( 18 ) 

i 

and consists of products of Wilson coefficients Cj and dimension-six operators Qj. The Wilson 
coefficients Cj and the matrix elements of the operators Qi both depend on the renormalization 
scale /i, which we set to the bottom-quark mass m b throughout this article when evaluating the 



5-meson observables of interest. In (18) the Fermi constant Gp and the leading CKM factor 
V t *Vtb have been extracted as a global prefactor. Finally, the sum over i comprises the current- 
current operators (i = 1,2), the QCD-penguin operators (i = 3,4,5,6), the electromagnetic 
and chromomagnetic dipole operators (i = 7, 8), and the semileptonic operators (i = 9, 10). 

For what concerns our work, the electromagnetic dipole and the vector-like semileptonic 
operator are the most important ones. We define them in the following way 

2 

Q7=j^m b (ri(sa^P R b)F flu , Q 9 =°(s 7 >>p Lb )(e 7(t e) , (19) 

where e is the electromagnetic coupling constant, m^/i) denotes the MS bottom-quark mass, 
is the field strength tensor associated to the photon, Pl,r — (1 T 7s)/2 project onto left- 
and right-handed chiral fields, and = i [7^,7"] /2. Within the SM one has numerically 
C7,SM(^b) ~ —0.3 and C 9i sm(^6) ~ 4.1. The chiral-flipped partners Q' 7 and Q' 9 of the electro- 



magnetic dipole and the vector-like semileptonic operator are obtained from (19) by simply 
replacing the chiral projectors Pl,r through Pr^l- 



In what follows, we supplement the effective SM Lagrangian ( 18 ) by 

i 

where A > m t denotes the scale of new physics and the index % runs over the complete set of 
dimension-six operators with flavor content (sb)(fr), namely (A,B = L,R) 

Qs,AB = (sP A b) (fP B r) , 

Qv,ab = (s^P A b) (f^P B r) , (21) 
Q T ,A = (s^ u P A b) (Tct^Pat) . 

As we will discuss in detail below, the Wilson coefficients of these operators can be bounded 
by various direct and indirect constraints for any given value of A. 

For later convenience, we also introduce Wilson coefficients C< associated to the non- 



standard operators (21) that are normalized like the SM contributions (18). In terms of the 



Wilson coefficients C± they are simply given by 

We now have enough definitions and notations in place to start the discussion of the relevant 
restrictions. 
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3 Direct Bounds on (sb)(rr) Operators 



The ten operators entering (20) govern the purely leptonic B s — > r 4 
X s t + t~ decay, and its exclusive counterpart B 



decay, 

K + T 



the inclusive 
making 



semi-leptonic B — > A s r 1 r aecay, ana its exclusive counterpart a 1 — > A 1 r 1 r 
these channels potentially powerful constraints. In practice, however, flavor-changing neutral 
current Bd, s decays into final states involving taus are experimentally still largely unexplored 
territory. Rather weak limits on the branching ratios of the former two modes can be derived 
from the LEP searches for B decays with large missing energy [38] and/or the 5-factory 
measurements of the ratio of the B s - and the I^-meson lifetimes [21]. In the latter case, one 
finds by comparing the SM prediction T B Jr Bd — 1 6 [—0.4,0]% [3] with the corresponding 
experimental result t Bs /tb<j — 1 = (—2.7 ± 1.5) % [7jJ^]the following bound 



B{B s ^t+t-) < 5%. 



(23) 



The same loose limit also holds in the case of BiB — > X s t + t ). Bounds of similar strength 
also follow from charm counting [391 HQ] . Combining the latest BaBar result for the average of 
the B- and B d charm multiplicity n c = 1.20±0.05 with B(B -> X scS ) = 0.24±0.02 @T], yields 
B(B — > X no c harm) = (4 ± 5)% at 68% CL. Subtracting from this number the SM prediction 

of around 2% 



for B{B — > X no charm) of around 2% (42], suggests that (23) is a conservative upper limit, 
that is consistent with charm counting^ This unsatisfactory experimental situation has been 
improved recently by the BaBar collaboration, which was able to set a first upper limit on the 
branching ratio of B + -)■ K + t + t~. At 90% CL this bound is 1441 



B(B' 



T T 



< 3.3-10- 



(24) 



Although the limits (23) and (24) are many orders of magnitude above the corresponding 



SM predictions, we will see below that they provide currently the strongest constraints on the 
Wilson coefficients of the scalar and vector operators of (s6)(fr) type. Let us first have a look 
at B s — > t + t~ and B — > X s t + t~ . In the case of the exclusive decay the relevant branching 
ratio reads jl5] 



B{B S 



— > T^T~ 



Ml 



1 - 



Am 2 T 



„ 2m r „ 



(25) 



where 



B s 



G 2 F M B J BjBa 
16^ lVM 



4 m 2 



6.5 • 10" 



(26) 



6 In fact, the current SM expectation and the measurement of the lifetime ratio differ by 1.5a [3]. If a 
significant deviation of Tg s / t g d — 1 from zero could be established, it would hint at new physics in T^ 2 . A 
more precise theoretical determination as well as an improved measurement of this quantity are thus important. 

7 The B — > X no charm branching ratio has also been determined by DELPHI 43 . The quoted 95% CL 
bound of B(B — > X no charm) < 3.7% relies however on a Monte Carlo simulation to model the displaced vertex 
distribution for intermediate charm, which introduces an uncertainty that is hard to quantify. 

8 The measurement is restricted to dilepton invariant masses of 14.23 GeV 2 < q 2 < (M B + — M K +) 2 . This 
experimental cut is implemented in our calculation. 
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and 

M Bs 



s,p 



m b + m s l 



Cs,RR T Cs,LL ± Cs,RL — Cs,LR 



01 

Fa = — ClO.SM + Cy,LR + Cy,RL — Cy,LL ~ Cy,RR 
Z7T 



(27) 



Here all Wilson coefficients are understood to be normalized as in (18) and (22) and evaluated 
at m b . Notice also that the SM contribution proportional to Ch^smC? 77 *) ~ —4.2 has been kept 
separately in Fa- Compared to the new-physics contribution this correction is suppressed 
by both a CKM and a loop factor, resulting in B(B S — > t + t~ )sm ~ 8 ■ 10~ 7 . In order to 
obtain this number, we have employed Gf = 1.16637 • 10~ 5 GeV -2 , M Ba = 5.3663 GeV [16"] . 
f Ba = 238.8 MeV g7], r Ba = 1.477 ps 0, \V*V tb \ 2 = 1.6 • 10" 3 @8], m T = 1.777 GeV, and 
a = a{M z ) = 1/129 [49]. 



Assuming the dominance of a single operator and neglecting the SM contribution to (27) 



we then find by combining (23) with (25) the following upper bounds on the magnitudes of 
the Wilson coefficients 

\CsM m b)\ < 0.7(0.6), \C VtAB (m b )\ < 1.4. (28) 

Here m b = m^ ole = 4.8 GeV (m b = m b {m b ) = 4.2 GeV) and m s = 0.1 GeV have been used. 
Notice that the choice of renormalization scheme employed for the bottom-quark mass is 
arbitrary given that we work at leading order (LO) in the perturbative expansion. We will 
use the differences of the results obtained for the two choices of m b given above to gauge the 
theory uncertainties plaguing our calculations. 

By means of the solutions to the LO renormalization group equations (RGEs) presented 



in Appendix |AJ the results (28) can be reinterpreted as limits on the matching corrections to 
the Wilson coefficients at the scale A. Performing the evolution from m b (m t ) up to m t (A) 
in a five-flavor (six-flavor) theory, we obtain in the case of the scalar operators 

|CWA)| = vl V23 vf \Cs,AB(m b )\ < V f 0.5 (0.4) , (29) 

where 775 = a s {m t )/a s {m b ) w 0.5 and r? 6 = a s (A)/a s (m t ) = [1.223 + 0.124 In (A/TeV)]" 1 . The 
bound on Cv,ab is, on the other hand, scale independent, because the vector operators Qv,ab 
correspond to conserved currents. In Table [TJ we summarize for convenience the currently 
available direct bounds on the high-scale Wilson coefficients assuming a new-physics scale of 
A = 1 TeV. 

Unlike B s — > t + t~ , the inclusive decay B — > X s t + t~ depends on all Wilson coefficients. 
Its differential branching ratio might be split into the SM contribution, the interference of 
the SM with the new-physics contributions, and the pure new-physics corrections. Since we 



assume that the Wilson coefficients of the operators (21) are generated at tree level, we are 
allowed to neglect both the SM and the interference terms, which are suppressed relative to the 
new-physics effects by at least a factor of order a/(4n) w 6 ■ 10~ 4 . The relevant contributions 



to the branching ratio differential in the dilepton invariant mass, s = q 2 , read [50 

I dB(B -> X.t+t-)\ 



ds 



') =A/ Mr+T - |a| 2 M^(s). (30) 

/ i 
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Ci(lTeV)\ 


B s -> T + T- 


B ->■ X s r+r- 


5+ -»■ K+r+r- 


S,AB 


0.4 


2.6 (5.8) 


0.5 (0.4) 


V,AB 


1.4 


2.1 (4.8) 


0.8 


T, A 




0.7 (1.6) 


0.5 



Table 1: Direct upper bounds on the high-scale Wilson coefficients at A = 1 TeV. The 
numbers without (in) brackets correspond to m b = 4.8 GeV (mj = 4.2 GeV). 



Here % = S, AB, V, AB, T, A and 



|V£K 6 | 2 £K£?^X c ^) cxp 



2m° l^csl 2 /(z)k(z) 



The analytic expressions for the kinematic functions Mj Xs (s) can be found in Appendix [B 
In order to obtain the numerical result in (31), we have employed | V^*Vt&| 2 /| V cb \ 2 = 0.96 |4"8~ 



B(B -> X c ^) exp = 10.23% 0, and z = (mf/mj*) 2 = 0.084. For this input the phase- 
space factor and the NLO QCD corrections jSI] of B — > XJLvi evaluate to f(z) ~ 0.54 and 
k(z) ~ 0.88. The Wilson coefficients are bounded by integrating the non-resonant branching 
ratio over the entire kinematical range 4m 2 . < s < (m b — m s ) 2 and comparing the obtained 
result to the experimental extractionn Numerically, the integrations yield 



dsM* s {s) = 2636 (170), 
J ds M* s (s) = 10542 (679) , (32) 

J dsM* a {s) = 126505 (8154) . 

From the bound B{B — > X s t + t~) < 5%, which is five orders of magnitude above the SM 
expectation B(B — > X s t + t~)sm ~ 5 • 10~ 7 [52], we then derive, by utilizing (30) to (32) and 
considering each type of Wilson coefficient individually, the inequalities 

\CsM^b)\ < 4.2 (9.6) , \C VtAB (m b )\ < 2.1 (4.8) , \C TA (m b )\ < 0.6 (1.4) . (33) 



Comparing these results with (28), one observes that B — > X s t + t~ constrains the contribu- 
tions arising from scalar and vector operators much less severely than B s — > t + t~ . 

It is again straightforward to find the corresponding bounds on the high-scale Wilson 



coefficients. In the case of the scalar operators, one proceeds in analogy to (29), the upper 
bound on Cv,ab is unchanged, and in the last case, we arrive at 

|CV,a(A)| = % 4/2 V /21 \C T , A (m b )\ < % - 4/21 0.7 (1.6) . (34) 

Experimentally the narrow ip 1 resonance at s « 13.7 GeV 2 has to be removed by making appropriate 
kinematic cuts in the invariant mass spectrum. In view of the poorness of the upper limit on B{B — > X s t + t~), 
it is at present immaterial if these cuts arc imposed in the theoretical calculation. 
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As in (31 ) to (33), the number in bracket corresponds to m b = 4.2 GeV, while the unbracketed 



value has been obtained for m b = 4.8 GeV. The bound (34) and the limits on Cs,ab and Cv,ab 
stemming from B — > X s t + t~ are collected in Table [I] 

Let us finally consider the bounds on the Wilson coefficients Cj that follow from the upper 



limit (|24|), which still leaves ample room for new physics to enhance the branching ratio with 

-7 



respect to the SM prediction of B{B + — > K + t + t~)sm ~ 2 • 10~' [53]. Keeping again only the 
pure new-physics contributions, the individual corrections to the differential branching ratio 
of B + -» K + t + t~ take the form [M] 



^dB{B + K+t+t- 



ds 



M, 



B+^K 



where 



Mi 



Gl\V*V tb \ 2 r B , 



B+^K+t+t- 



24vr 3 M| + 



+ T + T - ICA 2 M 



5.0-10" 



A' 4 



(35) 



(36) 



corresponding to t b + = 1.641 ps (7] and M B + = 5.279 GeV [15]. The functions M/ f+ (s) are 
collected in Appendix IbJ Integrating them over 14.23 GeV 2 < s < {M B+ - M K+ f gives 



dsMg + ( s ) 



dsM^ + {s) 



dsMfi + (s) 



Inserting (36) and (37) in to 
at a time, the upper bound (24) implies 

\C s ,AB(m b )\ < 0.8 (0.7) 



1082 (1422) 



1083. 



3614. 



(37) 



(35 ), and allowing for the presence of a single Wilson coefficient 



\C v ,AB{m b )\ < 0.8, \C TA {m b )\ < 0.4 



(3* 



T + T 



The limits on the Wilson coefficients at the scale A = 1 TeV following from B + — > K 
are shown in Table [Tj The quoted numbers have been obtained from (38) by simply apply- 
ing the relations (29) and (34). We see that the exclusive b — > st + t~ mode provides at 



present the most stringent direct constraints on the vector and tensor operators with flavor 
content (sb)(fr), while it is at pai 10 with B s — > t + t~ for what concerns the scalar contribu- 
tions. Notice also that the resummation of large leading logarithms of the form a" ln n (A 2 /ml) 
makes the bound on \Cs ; ab(A-)\ (|Ct,ab(A)|) relative to the limit that applies to \Cs,AB{™>b)\ 
(|Ct,ab(?^&)|) stronger (weaker). The running effects are however moderate in both cases if 
new physics enters at the TeV scale, as they change the results by a factor of around 1.6 and 
0.9 only. Numerically, we see that the direct constraints allow for effects that reach almost 
(9(1) in all |Cj(l TeV)|. Recalling that the dominant SM contribution to T\ 2 arises from the 
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One should keep in mind that the bound from B s 



is only an estimate, while the limit from 



B + — > K + t + t holds at 90% CL. In our numerical analysis of rj 2 , we will employ the most conservative 
bound on |Cs,ab(A)|. 
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color-singlet current-current operator (s7 a 'Plc)(c7 a1 Pi,6), which has a Wilson coefficient of 
roughly 1 at the weak scale, we expect that the (sb)(fr) operators can give a sizable correc- 
tion to the off-diagonal element of the decay matrix in the P s -meson system. We however 
curb our enthusiasm and postpone a detailed numerical analysis of the new-physics effects in 
2 to Section |5j to check first that the indirect constraints associated to operator mixing do 
not thwart these potentially large effects. 



4 Indirect Bounds on (s6)(rr) Operators 



Further constraints on the Wilson coefficients (|22j) arise indirectly from the experimentally 
available information on the b — > 57 and b — > s£ + £~ (£ = e, /i) transitions, because some of 



the effective operators introduced in (21) mix into 

2 

Qi,a = ~ 2 m T {sa^P A b)F, u , Q 9A = \ {sj»P A b)(£ % £) , (39) 
9s 9s 

at the one-loop level. Here A = L,R. The corresponding Feynman diagrams are depicted in 
Figure [ij It is important to realize that due to the flavor structure of the (sb)(fr) operators 
only insertions are possible in which the tau lines are joined and the photon is emitted from 
the resulting closed loop. Our results for the anomalous dimension matrix (ADM) describing 
the operator mixing of Qs,ab, Qv,ab-, and Qt,a into Q-j^a and Qg,A are given in Appendix |A~} 
We find that the ADM is sparse, because most of the relevant penguin diagrams either con- 
tain a vanishing Dirac trace or evaluate to zero due to current conservation. In particular, 
the operators Qs,ab mix neither into Q-j^a nor Qg^j^j] while Qv,ab (Qt,a) mixes only into 
Q9,a {Qi,a)- The possible one- loop mixing of the (s6)(rr) operators is therefore much more 
restricted than what has been claimed in the articles [27J [311 [33] , which all did not perform 
an explicit calculation of the one-loop ADM. We will see in a moment, that as a result of the 
particular mixing pattern, the stringent constraints from the radiative decay B — > X s j rule 
out large contributions to only if they arise from the tensor operators Qt,a- Similarly, the 
rare decays B — > X s £ + £~ and B — > K^£ + £~ primarily limit contributions stemming from the 
vector operators Qv,ab- 

The new-physics corrections to the partonic b — > sj and b — > s£ + £~ amplitudes involve the 



low-energy Wilson coefficients of the effective operators introduced in (39). Neglecting the 
matching corrections to both Qj t A and Qg t A, i-e., setting CV^A) = Cg^A) = 0, we find from 
the analytic solutions to the RGEs as given in Appendix |Aj the following expressions 

c 7 A^) = vt m (0.6 - % x ) C TjA (A) , 

(40) 

C 9A (m b ) = (0.1 - 0.2 r^ 1 ) (C VAL (A) + C v ,ar( A )) ■ 



11 Notice that scalar operators with flavor structure (sb)(qq) with q — s,b mix into the electromagnetic dipole 
and vector- like semileptonic operators at the one- loop level [55, 56. . The mixing arises from graphs constructed 
by joining two strange or bottom quarks belonging to the different disconnected parts of the composite operator 
and attaching the photon to the resulting open strange- or bottom-quark loop. 
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Figure 3: Diagrams with a penguin insertion of a (sb)(fr) operator (black squares) 
that contribute to the renormalization and the matrix element of the electromagnetic 
dipole (left) and vector-like semileptonic (right) operators. The tau loop in both graphs 
is closed. 



Determining the contributions of the operator sets (21) and (39) to the radiative and 
semileptonic b — > s processes also requires the knowledge of the corresponding tree-level and 
one-loop matrix elements. Details on these calculations are relegated to Appendix[Cj Following 
common practice, we include the effects of the relevant matrix elements by defining so-called 



effective Wilson coefficients [57]. For the SM operators Q 7 and Qg appearing in (19), we obtain 
at the scale \i = mj, the following new-physics corrections 

ACf(s, m b ) = y/x~ T (C 7tR (m b ) + M 7 (s, x T ) C T)R (m b )) , 

(41) 

ACf(s,m b ) = C 9 , L (m b ) + M 9 (s,x T ) {C v , LL (m b ) + C ViLR (m b )) , 



with C 7) A( m b) and C 9t A(m b ) given in (40) and CV,A( m &) = 0.9 r]^ 21 Ct,a{&). Furthermore, 
s = q 2 /m 2 and x T = ml/m 2 . The analytic expressions for the functions M 7 (s,x T ) and 
Mg(s,x T ) can be found in Appendix [C] Analogous formulas hold in the case of the Wilson 
coefficients AC 7 cS (s,m b ) and AC' 9 cG (s, m b ) of the chiral-flipped operators Q' 7 and Q' 9 after 
replacing Ci^{m b ) through C^L(m b ) (Ci^A^b) through C^RA{. m b)) in the first (second) line 



of {41). Notice that only ACf(m b ) = ACf(0,m b ) and AC' 7 eS (m b ) = AC' 7 cS (0,m b ) enter the 
prediction for B — > X s 7, while all AC^ s (s,m b ) and their chirality-flipped partners affect the 
b — > s£ + £~ transitions. 

In order to derive a bound on the Wilson coefficients Ct,a at the high scale, we compare 
the experimental measurement of the B — > A s 7 branching ratio with its SM prediction. For 
a photon-energy cut of E n > 1.6 GeV, the experimental world average reads [7J 

B(B -»■ A s7 ) oxp = (3.55 ± 0.24 ± 0.09) ■ 10" 4 , (42) 

while the SM expectation is given by [581 ESPf] 

B{B X s7 ) S m = (3.15 ± 0.23) ■ 10" 4 . (43) 

"Several NNLO QCD corrections (see [60l EU |62l EH ISU |65] and partly EBJ) that were calculated after 
the publication of 58, 5H] are not included in the central value of the SM prediction, but remain within the 
perturbative higher-order uncertainty of 3% that has been estimated in the latter two articles. 
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The good agreement between (42) and (43) puts stringent limits on the corrections ACf ff (m&) 



and AC' 7 cS (mb). Further non- negligible constraints on the magnitudes and phases of the 
corrections to the low-energy Wilson coefficients of the electromagnetic dipole operators arise 
from the inclusive [67] and exclusive [68] b — > s£ + £~ transitions as well as the time-dependent 
CP-asymmetries and the isospin asymmetry in B — > K*j |69j . Including all these constraints 
and allowing for the variation of a single effective Wilson coefficient at a time, we obtain the 
following 90% CL upper bounds 

|AC 7 cfr (m 6 )| < 0.29, |AC; cff (m 6 )| < 0.19. (44) 

Notice that since we are treating the Wilson coefficients as complex, the upper limit on the 
magnitude of ACj R (mb) is weaker by almost a factor of 6 than the bound that holds in the 
case of a real coefficient. In contrast, the bound on \AC' 7 eS (mb)\ is not affected by whether the 
correction is treated as complex or real, because Q' r does not interfere with the SM contribution 
to first order. We also remark that the B — > K^£ + £^ observables are more efficient than 
B — > K*^ in restricting the allowed values of the real and especially the imaginary parts of 
ACf" (///,,). 



Utilizing (40) and (41), the limits (44) translate into 



\C T>R (A)\ = rff' 21 (0.4 + 2.3 ^y 1 \ACf(m b )\ < % 4/21 (8.1 + 1.4% T' , (45) 

and an analog inequality, obtained by replacing \ACj S (rrib)\, 8.1, and 1.4 by \AC' 7 eS -{mb)\, 12.4, 
and 2.2, holds in the case of Ctl- F° r a new-physics scale of A = 1 TeV the numerical values 
of the bounds on |CV_a(A)| are given in Table |2j Notice that as a result of the resummation 
of large logarithms the bounds on the low-scale Wilson coefficients of the tensor operators are 
a factor of around 1.2 stronger than those on the initial conditions. 

The magnitudes of the Wilson coefficients CV,ab can be bounded by comparing the available 
data on the inclusive [701 EJ E21 [73] and exclusive [Hll5l[76l[77J[79l[78l|80]6^ s£ + £~ transi- 
tions with the corresponding theoretical predictions that include the corrections ACg S (s,rrib) 
and ACg (s, nib). Our global analysis of rare semileptonic B decays relies on the article [1H] for 
what concerns B X s £ + £~ and on the works [681 EH E21 E3J EI] in the case of B -> K^£ + £~. 
The used dataQ includes the branching ratios B(B -»■ X s £ + £~) and B{B K i *H + £-), the 
forward-backward asymmetry Afb{B — > K*£ + £~), as well as the longitudinal K* polarization 
fraction F^(B — > K*£ + £~). In all cases, we consider the low-s region of dilepton invariant 
masses s = [1, 6] GeV 2 . At high s we consider for B — » X s £ + £~ the region s > 14.4 GeV 2 and 
for B ->■ K { *k + £~ the two bins s = [14.18, 16] GeV 2 and s > 16 GeV 2 . At 90% CL our global 
analysis returns the following upper bounds 

\ACf(s,m b )\ < 2.0, \AC' 9 c{i (s,m b )\ < 1.7, (46) 



13 Our fit is based on [7Ql [7T] , [22 EH [79l [80] , and [76l [78] for what concerns B -> X s l + £~, B -> K*l + l~, 
and B — > K£ + £~ , respectively. The latest (but still preliminary) Belle measurement of B(B —> X s £ + £~) [72] . 
which finds a low-s branching ratio of around 2.4<r below the SM prediction, is not used, because the quoted 
central values of the branching fractions for electrons and muons in the final state differ by a factor of more 
than 2. If only the LHCb measurements of B -» K*£ + £~ [79] are used, the numbers in (44), (46), and Tabled 



read 0.37 and 0.19, 3.1 and 2.5, and 1.5, 1.2, 0.07, and 0.16, respectively. While some of these bounds are not 
as strong as the limits based on the combined B — > K^*H + £" data [751 ESI 1ZH1 [HO] , they clearly show the large 
impact (and future potential) of LHCb measurements on the extraction of low-energy Wilson coefficients. 
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Figure 4: 1PI diagram with a penguin insertion of a (sb)(fr) operator (black squares) 
that contribute to the b — > 377 amplitude. The tau loop in the graph is closed and the 
diagram with interchanged photons is not shown. 



when the corrections to the effective Wilson coefficients are treated as additive, s-independent 
contributions and only a single correction is allowed to float at a time. 



The inequalities (46) can be converted into bounds on the magnitudes of the high-scale 



values of Cv,ab using (40) and (41). Since the expressions for ACg S (s,mb) and AC 9 eff (s,mb) 
depend explicitly on s through the function M 9 (s, m T ), we extract the limits directly from our 
global fit. We find the following exclusion 

\CvM A )\ = ( 2 -° " O- 2 %T'l AC f(s,m b )\ < 2.0 (2.0 - 0.2 r^ 1 )" 1 . (47) 

An analog formula holds for |CV,_ra(A)| after replacing |AC| ff (s, m.&)| and the factor 2.0 in 
front of the bracket by \ACg cS (s, m&)| and 1.7. Assuming a new-physics scale of A = 1 TeV, 
one finds the numerical values given Table [2} 

In contrast to B — > X s 7, all (s6)(fr) operators contribute to the double-radiative B ds - 
meson decays at the one- loop level [Sol 186"] . In the following, we will concentrate on the case 
of B s — > 77, which turns out to give the strongest constraints on the Wilson coefficients of the 
operators under consideration. At the quark level the double-radiative b — > 577 decays receive 
contributions from one-particle irreducible (1PI) diagrams of the type shown in Figure [4] as 
well as similar one-particle reducible (1PR) graphs in which at least one of the two photons is 
radiated from an external leg. Since the 1PR diagrams are proportional to the tree-level matrix 
elements of Q?\ their contributions are included if the corresponding amplitudes are expressed 
through the effective Wilson coefficients C^ eS . A detailed discussion of the calculation of the 
one-loop matrix elements can be found in Appendix [D| 

Experimentally so far only upper limits on the branching ratio of B s — > 77 exist. At 
% CL the most stringent limit is [HZj 



B{B S 77 ) cxp < 8.7 ■ 10~ 6 . (48) 

In order to constrain the Wilson coefficients of the (s6)(rr) operators also a value for the 
branching ratio of B s — > 77 within the SM is needed. Updating the analysis of [881 EH] , we 
obtain 

B(B s ^ 11 ) SM =(0.7tH)-lO-\ (49) 



15 



|Q(lTeV)| 


b — > sj, s£ + £ , S77 


S,AB 


— 


V,LA 


1.1 


V,RA 


1.0 


T,L 


0.07 


T,R 


0.10 



Table 2: Indirect upper bounds on the high-scale Wilson coefficients at A = 1 TeV. 



where the dominant part of the quoted error is due to the hadronic parameter A#, which 
parametrizes the first negative moment of the B- meson wave function. Estimates of A# are 
very uncertain, but typically fall in the range between 0.25 GeV and 0.75 GeV [9"U t I9~T | I9"2"| I9"5] . 



and the central value in (49) corresponds to \b = 0.5 GeV. Numerically subleading errors in 
the latter SM prediction arise from the variation of the renormalization scale and the error on 
the decay constant fs s . 



Assuming the dominance of a single operator, we then find by combining (48) and (49) the 



following 90% CL upper bounds on the magnitudes of the Wilson coefficients (A, B = L,R) 

\ACf{m b )\ < 2.2, \AC' 7 eS (m b )\ < 1.9, 

(50) 

\C SAL {m h )\ < 3.4, \C s>AR {m b )\ < 2.3, \C V)AB {m b )\ < 5.9. 

Clearly, these limits are not competitive with the bounds obtained earlier from the other tree- 
and loop-level mediated 5rf iS -meson decays. The reason for the weakness of these bounds is 
threefold. First, as shown in Appendix IdI the 1PI contributions of the (sb)(fr) operators are 



power suppressed, second, the experimental bound (48) leaves room for order of magnitude 
enhancements of B(B S — > 77) with respect to the SM expectation, and third the theoretical 
errors plaguing ( |49~| are large. Notice that the power suppression of four-fermion operator 
contributions is a generic feature of the B s — > 77 decay in the heavy-quark limit [HSl I8"9"] . 
which is absent in the case of B — > ^77 [86J. The sensitivity of the latter decay mode to 
(s6)(tt) operators is therefore in principle more pronounced than the one of B s — > 77. In 
practice, however, the possibility to extract short-distance information from B — > ^77 is 
severely limited. First, the perturbative part of the B — > K'-yj spectrum is not accessible 
behind the large resonance peaks associated to B — > Kr\ c — > K^yy etc. and, second, not even 



an upper limit on B[B — > K'yy) is currently available. Hence the bounds (50) represent at 
present the most stringent limits on the low-energy Wilson coefficients of interest that can be 
derived from double-radiative -B^-meson decays. 

With a large amount of luminosity collected at a super-flavor factory it might also be 
possible to measure the CP asymmetry 



_ l4(^77)| 2 -lA(i^77)| 2 
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besides the B s — > 77 branching ratio. The SM prediction of rep reads 

(r C p)sM = (0.5+°J)%. (52) 
The largest fraction in the quoted error is due to the scale dependence, followed by the uncer- 



tainty coming from A#. The smallness of (52) implies that finding tq P at or above the 10% 



level would constitute an unambiguous signal of the presence of additional sources of direct CP 
violation. In fact, the complex Wilson coefficients of the (sb)(fr) operators give rise to such 
an effect and thus allow to change vq P dramatically. Employing the 90% CL bounds on the 
low-energy Wilson coefficients of \C StAB {m b )\ < 0.8, \C VAB (m b )\ < 0.8, \ACf(m b )\ < 0.29, 
and \AC' 7 eS (mb)\ < 0.19, we find that the CP asymmetry in B s — > 77 can vary in the 
ranges (r CP ) S:AL G [-95,95]%, (r CP )s,AR e [-40,40]%, (r CP ) v>AB G [-70,70]%, (r CP ) T ,L G 
[—20, 20]%, and (r CP ) TjR G [—70, 70]%, for the respective Wilson coefficient with an arbitrary 
phase. In order to measure r CP , SuperB (Belle II) probably needs to collect at least 75 ab -1 
(50 ab -1 ) of data, a goal that is expected to be achieved in the year 2021. This implies that 
the B s — > 77 decay will not be able to provide useful additional information on the possible 
size of (s6)(fr) contributions in the near future. 

By comparing the results displayed in Table [T] with those shown in Table [2| one observes 
that the bounds on the high-scale Wilson coefficients Cj(A) that stem from the direct con- 
straints, i.e., the rare decays B s — > t + t~ , B — > X s t + t~ , and B + — > K + t + t~ , are currently in 
8 out of 10 cases stronger than those that follow from the indirect constraints, i.e., the b — > 37, 
b — > s£ + £~ , and b — > 577 transitions. The only exception are the Wilson coefficients Ct,a(A) 
of the tensor operators, which are tightly bound as a result of the one-loop mixing of Qt,a into 
the electromagnetic dipole operators Q-j and Q' 7 , which give the dominant contribution to the 
B — > X s 7 rate. On the other hand, sizable effects in Cs^^A) and Cv,ab(A), that can almost 
compete with the leading current-current SM interactions in strength, are allowed, because 
rare B s and Bd decays involving a tau pair in the final state are so poorly constrained. While 
the direct constraints are unlikely to change notable with almost all BaBar and Belle data 
analyzed, the indirect bounds on the vector operators are expected to improve with LHCb 
precision measurements of B — > K^*H + £~ under way. 



5 Effects of (s6)(tt) Operators in 

The off-diagonal element of the decay-width matrix is related via the optical theorem to the 
absorptive part of the forward-scattering amplitude which converts a B s into a B s meson. 
Working to LO in the strong coupling constant and Aqcd/^6, the contributions from the 



complete set of operators (21) to T\ 2 is found by computing the matrix elements between 
quark states depicted in Figure |5j Details on the actual calculation are given in Appendix [Ej 
Considering again only the self-interference of a single (sb) (f r) operator at a time, we obtain 
in terms of the hadronic matrix elements (Qi) = (B s \Qi\B s ) of the following AB = 2 operators 

Qj = (sP A b) (sP A b) , 

Q A s = {s a P A b p ){spP A b a ) , (53) 
Q$ = (8-fPAb) {si,P A b) , 
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Figure 5: Diagram with an double insertion of a (sb)(fr) operator (black squares) that 
contributes to rf 2 at LO. The tau loop in the diagram is closed and the cut (dashed 
line) indicates that only the imaginary part of the graph furnishes a correction to the 
off-diagonal element of the decay-width matrix. 



the results 



(r 



12) S, LA 



Af r s 2 3x T f3 T (Q L s )(C s ,LA) 2 , 



x T ) (Q*) + (1 + 2x T ) (Qf > (C ViLA ) 



fil2JV,LA ~ 

(r? 2 ) T)£ = AT r?2 12x r /3 r [ (Q L S ) + 2 (Q L S ) j (C t ,l) 2 • 
Here a, f3 are color indices and we have introduced (3 T = a/1 — 4x T as well as 



(54) 



n 2 m 2 

F b (K s v tb y 



6vrM 



(55) 



The corresponding expressions for (rf 2 ) s RA , (T^vrai an< ^ 0^12)tr are obtained from (54) 
by exchanging the labels L and R everywhere. Following common practice, we express the 
matrix elements in terms of hadronic parameters Bi. In our analysis, we employ 



(Qs) 
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f B MlBs , (Qj) = y 2 flMlB s , (Q$) = \ f B Ml s B, 



(56) 



with f Bs = 231 MeV 0, B s _« 1.3, B s « 1.4 |SS], 5 y « 0.84 [S], and the meson states 
normalized as (B S \B S ) = (B S \B S ) = 2Mb s . Notice that the latter value of fs s has been used 
to obtain the SM predictions (Ar s )sM and (g^Jsm presented in rtij) and (J7J), respectively. 

The direct and indirect constraints discussed in Sections [3] and [4] restrict the magnitude 
of the low-energy Wilson coefficients but not their relative weak phase. We can thus deter- 
mine only the maxima^] allowed value of the parameter i?p introduced in ( 14 ) in a model- 
independent fashion. Allowing the new-physics contribution of one operator at a time to be 



14 Of course, there is also a model-independent lower limit on (i?r)i- It corresponds to full destructive 
interference and is obtained from (57) and (58 1 by replacing the plus with minus signs. 
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Figure 6: Left (Right): Constraints from Ar s and a s j- s on the parameter R? and 0r 
in scenario S2 employing the data set Dl (D2). The yellow and gray bands represent 
68% CL regions (dofs = 1), while the purple ellipses illustrate the parameter space that 
is accessible for \Cs,AB{jnb)\ < 0.8 and Cv,ab{th^) = Ct,a(?™&) = 0. For comparison, 
the 68% (95%) probability region of the combined fit to the data (dofs = 2) is shown 
in red (light red). In both panels the SM prediction (best-fit solution) is indicated by 
a dot (cross). 



non-zero, the corresponding inequality is given in terms of the off-diagonal elements of the 
decay matrix (54) by 

(Rr)i < 1 + cos (^)sm « 1 + . (57 ) 



Numerically, we find for the individual contributions 

(Rt)s,ab < 1 + (0.4 ± 0.1) \CsM m b)\ 2 , 

(Rf)v,ab < 1 + (0-4 ± 0.1) \C v M m b)\ 2 , (58) 

(R r ) TA < 1 + (0.9 ± 0.2) \C TA (m b )\ 2 , 

where the quoted uncertainties are due to the error on (Ar s )sM as given in Q. Employing 
now the 90% CL bounds on the low-energy Wilson coefficients derived in the previous two 
sections, i.e., \Cs,AB( m b)\ < 0.8, \Cy,AB( m b)\ < 0.8, |Ct,l(^&)| < 0.06, and |CV,B( m ;>)| < 0.09, 
it follows that 

(Rt)s,ab < 1-4 , (Rt)v,ab < 1-3 , (Rt)t,l < 1-004 , (R T ) T>R < 1.008 . (59) 
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These numbers imply that (sb)(rr) operators of scalar (vector) type can lead to enhancements 
of I rf 2 1 over its SM value by 40% (30%) without violating any existing constraint. In contrast, 
contributions from tensor operators can alter |rf 2 | by at most 0.8%, because larger corrections 
would be at variance with B — > X^. 

In Figure [6] we compare the parameter space in the Rr cos (pr^Rr sin 0p plane that can 
be populated in scenario S2 with \Cs,AB{mb)\ < 0.8 and Cv,ABijnb) = Cr,A(rnb) = (purple 
ellipses) to the constraints imposed by the measurement of Ar s and a s j s (yellow and gray 
bands). The left and right panel correspond to the data set Dl and D2, respectively. We 
see that the contribution of a single scalar operator with flavor structure (sb) (rr) leads to an 
improvement of the fit. For example, in the case of the new data set D2 an agreement with the 
data at 68% CL (red area) can be achieved. It is however also evident that the best-fit solution 
(cross) of both data sets cannot be accommodated. The same statements also apply to the 
case of vector operators, which we do not explicitly display in the figures as the accessible 
parameter space essentially resembles the one shown. This demonstrates that absorptive new 
physics in in form of (sfe)(rr) operators is only able to partly resolve the anomaly in the 
dimuon charge asymmetry data observed by the D0 collaboration. 



6 New-Physics Models 

So far our discussion has been completely general since we parametrized the new-physics effects 
in terms of effective couplings of higher-dimensional operators. We now will consider explicit 
scenarios of new physics that can give rise to (sb)(fr) operators with large Wilson coefficients. 
Discussions similar to ours have been presented previously in J2U [301 EH ES] . 

We first consider scenarios with leptoquarks (LQ), i.e., color-triplet or -antitriplet par- 
ticles carrying both lepton (L) and baryon (B) number. Such new degrees of freedom are 
encountered in various extensions of the SM, such as grand unified theories, technicolor, and 
composite models (see (961 E3> EHl ES] for topical reviews). At low energies, the interactions 
between the SM particles and the LQs can be captured by writing down the most general 
Lagrangian compatible with renormalizability, L and B conservation, and invariance under 
the SM gauge group. These general requirements allow for scalar and vector LQs interacting 
directly with the SM fermions as well as for LQ couplings to the Higgs boson. 

In the following, we will elaborate in detail only on the case of 577(2) singlet scalar LQs. 
From our discussion it should however become clear how the given results have to be adapted to 
cover the other cases. The relevant interaction terms involving a charged lepton, a down-type 
quark, and a 577(2) singlet scalar LQ are given by 

£ L q D (\ R s )ij d C jP R ei S + h.c. . (60) 

Here A R ^ o is a complex 3x3 matrix in the lepton and quark flavor spaces, di and ej are SU(2) 
SM singlets, and the subscript of So indicates the £77(2) transformation property of the LQ. 



After integrating out the LQ and performing a Fierz rearrangement, the interactions (60) give 
rise to a AB = 1 (sb)(fr) operator of vector type. Explicitly, one finds 

^cff 3 ^2 (°cV,RR j (OlJ 

So 
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where M§ denotes the mass of Sq. The coefficient multiplying the operator Qv,rr is bound 
by the requirement \Cv,ab{h)\ < 0.8, which holds for any scale \i. Numerically, we find 



M l 



< 2.1 TeV" 



(62) 



The same bound also applies to the possible SU(2) doublet scalar LQs (commonly denoted by 
62 and S2 m the literature), which generate Qv,lr operators, after an appropriate replacement 
of couplings and masses. 

Importantly, the interactions (60) give rise to corrections to both rf 2 and Mf 2 . The 
contribution to the former quantity arise from a double insertion of (61 ). Applying the general 
specific case of a sc 

(( A W 32 ( A ^ ) 33 )' 



formulas (54) to the specific case of a scalar LQ, we obtain 

, 2 



(r 



12)LQ 



288tt 



M- 



0r 



l-x T )B v --(l+2x T )Bc 



mif B M Bs . (63) 



This formula disagrees with the analytic expression first given in [36J , which does not contain 
terms proportional to B$ at all. Numerically, the difference amounts to a factor of around 
—0.4, with our result being smaller in magnitude. The correction to M{ 2 originates from the 
box diagram displayed on the left-hand side of Figure [7j In agreement with [SB] , we find for 
this contributionE3 



(M: 



12/LQ 



fj B J(a T )f 2 B M Ba B, 



(64) 



where fj Bs = %^ 23 p= 
function 1(a) reads 



384vr 2 

0.8 encodes the QCD corrections due to RG running and the Inami-Lim 



1(a) 



1 + a 



+ 



2a 



1-a) 3 



In a . 



(65) 



Notice that the above function approaches 1 for a — > 0, which is the relevant limit, since 
a T = mtlM\ 1 for LQ masses at or above the electroweak scale. 



Combining now (63) and (|64h , one obtains numerically 



12JLQ 



2084 



Ms 



six correction (sing 



250 GeV 

(^12) sm/ (rf 2 )sM 



2084 z 



So ' 



(r? 2 )LQ 

which is real and positive. In contrast, r§M 

Notice also that (66 ) scales as M| which reflects the fact that (Mf 2 )LQ arises from a dimension 



-200 exp (i(f> s su ) 



(66) 
-200. 



e operator insertion), whereas (r^ 2 )LQ stems from a dimension-eight con- 
tribution (operator double insertion). 

In fact, given that AM S (or equivalently Rm) is bounded experimentally, the sign and 
magnitude of tlq completely determine the ranges of possible Rr, AT S , and a s j s values in the 

15 We integrate out the top quark and the LQ simultaneously and consequently neglect small RG effects in 
the six-flavor theory. 
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b Sq s 




Figure 7: Left: One-loop contribution to Mf 2 arising from a box diagram involving tau 
leptons and SU(2) singlet scalar LQs. Right: Tree-level contribution to M{ 2 associated 
to the exchange of a Z' boson. 



LQ scenario (60). In Appendix [Fj we show that this is a feature of all new-physics scenarios 
with real (M^ 2 ) NP /(r^ 2 ) NP . Applying the general formula (F4) to the case at hand, we find at 
90% CL 

n n? n 9 

0.96 = 1 - — < (R r ) LQ < 1 + — = 1.31 . (67) 



'So 



So 



Here the bound on Rm given in (16) and (17) has been used and the final numerical values 



have been obtained for M$ o = 210 GeV, which represents the current 95% CL lower bound on 



third-generation scalar LQs decaying to a br final state |100j . Comparing (67) to the relevant 
model-independent upper bound, i.e., the inequality in (59) involving (R?)v,ab, we see that at 
present both limits are (accidentally) almost the same. Since the bound on M§ will improve 
with the LHC collecting more and more luminosity as we speak, the model-dependent upper 
limit will, however, soon become stronger than the model-independent one. 



Similarly, one finds from the general formulas (F5) and (F6), the following two double 
inequalities 



-0.11 ps^ 1 



-0.7-10" 



2\ / 02' 

if — ] (Ar s ) S M < (Ar s ) LQ < i + — ) (Ar s ) SM = o.ogps- 1 



z 



So 



z 



S 



(68) 



-331 1 + 



0.1 



Z So 



i/JsM < KJlq < 331 1 



0.1 



Z S 



0.7- 10" 2 . 



These numbers imply that light singlet scalar LQs can change the predictions for Ar s and 
a S f s with respect to the SM expectations Q and (j7j) by a factor of more than —1.3 and ±375, 
respectively. In contrast, a strong enhancement of Ar s is not possible. The maximal effects do, 
however, not occur simultaneously, since the observables in question are strongly correlated, 
as can be seen from Figure [8} In both panels the parameter space that is accessible at 90% CL, 
assuming a singlet scalar LQ mass of 210 GeV (400 GeV), is indicated in yellow (green). We 
see that Ar s becomes maximal (minimal) for 0}/^ = 0° (0jm<£ = —180°) and that the LQ 
predictions form a band with a cosine-like shape (left panel). The predictions for a^ s , on the 
other hand, are minimal for (j> s ji^ = 0°, —180° and maximal for 0}/^ = —90°, and the borders 
of the accessible parameter space in the 4>j/M~ aS fs plane are sine curves (right panel). Notice 



that the cosine- and sine-like behavior of the predictions is an immediate consequence of (15). 
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-150 -100 -50 50 



Figure 8: Predictions for Ar s (left) and a s j s (right) as a function of 0}/^ in the presence 
of a singlet scalar LQ with M§ = 210 GeV (yellow areas) and M§ o = 400 GeV (green 
areas). Since the left (right) panel is symmetric (antisymmetric) under (fij/^ -H- —^j/^ 
only positive values of 4> s j/^ up to 50° are shown. In both panels, the blue (red) error 
bar corresponds to the data set Dl (D2), while the SM point is indicated as a black 
point. See text for details. 



By comparing the accessible parameter space to the experimentally preferred regions, we see 
from the yellow colored regions that at 68% CL both the data set Dl (blue error bars) as well 
as D2 (red error bars) cannot be accommodated by a light singlet scalar LQ. Increasing the 
LQ mass, reduces the allowed parameter space, as can be seen from the green areas, further 
limiting the possible improvement of the tension in the 5 s -meson sector, 
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The above general considerations are readily applied to the case of vector-like LQs. Consid- 
ering for example LQ interactions of the form £lq D (^RV )tj djl^PR^-i Vj^ + h.c. will lead to an 
effective Lagrangian £ e ff 3 — {^rv^sz^rvo)^/ My o Qv,rr with My denoting the mass of the 
vector-like LQ. It follows that the model- independent bound on such a LQ is twice as strong as 



the limit derived in (62). Furthermore, the ratio between (M^ 2 ) L q and (r^ 2 ) L q is found to be 



equal to the result given in (66). This tells us that (for M Vo = M§ o ) the possible effects in Ar s 
and a s j- s due to V are identical to those originating from S . Vector-like LQs thus fail to repro- 
duce the data sets Dl and D2 at 68% CL as well. Notice finally, that the situation is utterly 



16 While this work was being completed, the preprint [101J appeared, which performs a dedicated global fit 
of scalar LQ couplings, utilizing an assortment of tree- and loop-level observables in the charged lepton and 
quark sector. This article finds that the tension in the B s -meson system cannot be cured by LQs, if one wants 
to explain the anomaly of the anomalous magnetic moment of the muon, (g — 2)^, and simultaneously fulfill 
the bounds on lepton-flavor violating processes. Notice, however, that if the former requirement is dropped, 
no bound on \{^* R § ^(A^g^sl/M- can be derived from the rare decays considered in the latter paper. 
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hopeless, if one considers LQ couplings to the Higgs boson, since after a Fierz rearrangement 
the corresponding effective scalar interactions with flavor structure (sr)(rfe) contain tensor 



operators, which are most tightly constrained by the model- independent bounds (59) 



We now switch gear and consider SM extensions where the gauge group contains an ad- 
ditional U(l)' factor and the resulting Z' boson possess family non-universal, flavor-changing 
couplings [102] . Such new degrees of freedom arise in many well-motivated new-physics mod- 
els, such as grand unified theories and/or string constructions (see [4"6 l I103j for comprehensive 
reviews). In the physical basis, the flavor-changing neutral currents generically appear at tree 
level in both the left- and right-handed sectors. The interaction Lagrangian relevant for the 
further discussion can be written as 



9 



cos 9 



(^ b sYPLb + K? b sYPnb + h.c.) + 



,L - 



k~ t t 7- P L t + /<; r r 



Prt 



Z' (69) 



where g is the SU(2) coupling of the SM and cos 8w denotes the cosine of the weak mixing 



angle. The Z'-boson coupling constant does not appear explicitly in (69), because it has been 



absorbed into the k-' factors. Notice that the flavor off-diagonal couplings are in general 
complex, while the diagonal ones have to be real due to the hermiticity of the Lagrangian. 
Integrating out the Z' boson leads to the following AB = 1 four-fermion interactions 



£ C fT D 



8G F Mf A B 



V2 M 2 zf 



K sb K TT Qv,AB 



(70) 



Here Mz> denotes the mass of the extra U(l) ; gauge boson. From the numbers collected in 
Table [TJ it follows that 

, . ) , , Li 

(71) 



I I 

1 sb rrl < 1.9 TeV" 



The interactions (69) give rise to tree- level corrections to both r* 2 and Mf 2 . In the fol- 



lowing, we restrict ourselves to the case of purely left-handed interactions. The inclusion of 
right-handed currents is however straightforward. For r^ 2 , we find in accordance with [30] 



(T 



vi)z> 



4G 2 F ( M z 



2 \ 2 



9tt \Ml, 



{ K sb K Tr) Pt 



[l-x T )B v --(l + 2x T )B s 



m 2 J 2 B M Bs . (72) 



The element M{ 2 receives contributions from the graph shown on the right-hand side of Fig- 
ure [7j A simple tree- level calculation leads to |104j 



(M? 2 )z' = ^^K) 2 VbJIM b B 



v ■ 



(73) 



The ratio of ( 73 ) and ( 72 ) is given in semi- numerical form by 

2 



r z > = 6.0 • 10 J 



M 7 , 



250 GeV kL 



6.0 • 10 5 y z , . 



(74) 



It is real and positive. Notice that the numerical coefficient in (74) is enhanced relative to 
(66) by a loop factor and thus very large. From the discussion in Appendix [F] it hence follows 
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Figure 9: Predictions for Ar s (left) and a s j s (right) as a function of 0}/^ in the presence 
of a Z' with M Z ' = 400 GeV (yellow areas) and M z i = 1000 GeV (green areas). Since 



the left (right) panel is symmetric (antisymmetric) under 



only positive 



values of 4>j/^ up to 50° are shown. In both panels, the blue (red) error bar corresponds 
to the data set Dl (D2), while the SM point is indicated as a black point. See text for 
details. 



that the effects of a Z' boson in R^, Ar s , and a s j s are severely constrained. Explicitly, we find 

1 1.0 ■ 10- 4 ,„ , 1 7.7- 10- 4 
1 < (Rr)z> < 1 + , 

yz< yz> 

( 7.7-10~ 4 \ / 7.0-10~ 5 \ 

- i + (Ar s ) SM < Ar s ) z , < i + (Ar,) S M , (75) 

V vz> ) \ yz> J 

f SS-10 _4 \ / 3S-10~ 4 \ 
-331 1 + (oJJsm < ifi).) z , < 331 1 + (aJ.)sM • 

V yz> J V yz> J 

The Z'-boson mass M z > is constrained by direct searches for resonant production of tau 
lepton pairs. At present the best bound is provided by the CDF measurement |105j . which 
rules out a Z' with SM-like couplings to qq and t + t~ and a mass below 399 GeV at 95% CL. 
Although this bound is model-dependent Q we will use it as a guideline and consider a rel- 
atively light Z'-boson of M z > = 400 GeV as well as a heavier one with M z > = 1000 GeV. 
The Z'-boson predictions in the 0j,^-Ar s ((fij/^-a'fs) plane are shown in the left (right) 



17 The mass Mz 1 and the couplings k^ r are constrained by the LEP measurement of the Z-boson couplings 
to taus [106] . Applying the results of the one-loop Z —> ff form factors given in |107) , we find that even for 
couplings close to the non-perturbative limit the resulting bound is however with Mz 1 ^ 25 GeV rather weak. 
We will not entertain the possibility of a very light Z'-boson in the following. 
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panel of Figure |9j The parameter space of 90% probability is depicted in yellow (green) for 
Mz> = 400 GeV (Mz> = 1000 GeV). While the shapes of the colored areas resemble the ones 
of Figure |8j one observes that compared to the case of LQs the allowed parameter space is 
even more reduced for a U(l)' gauge boson. In particular, a Z' boson with mass at or above 
the electroweak scale and purely left-handed perturbative couplings clearly fails to describe 
the 5 s -meson data within 68% CL. Because the parameter rz> as defined in (74) is of the 
order m^ 2 M\, j M\ and therefore generically large, the latter statement is also true if the 
Z' boson couples only right-handed or both left- and right-handed. 



7 Summary 

Motivated by the observation that the anomalously large dimuon charge asymmetry measured 
by the D0 collaboration, can be fully explained only if new physics contributes to the absorp- 
tive part of the B s -B s mixing amplitude, we have presented a model-independent study of 
the contributions to arising from the complete set of dimension-six operators with flavor 
content (sb)(fr). Taking into account the direct bounds from B s — > t + t~, B — > X s t + t~, 
and B + — > K + t + t~ as well as the indirect constraints from b — > 57, b — > sl + l~ (£ = e, /x), 
and b — > 377, we have demonstrated that only the Wilson coefficients of the tensor operators 
are severely constrained by data, while those of the scalar and vector operators can be sizable 
and almost reach the size of the Wilson coefficient of the leading current-current operator in 
the standard model. The model-independent 90% CL limits on the magnitudes of the Wilson 
coefficients are summarized in Table [3j Given these loose bounds, it follows that the presence 
of a single (s6)(fr) operator can lead to an enhancement of rf 2 of up to 40% compared to 
its standard model value, thereby reducing but not entirely resolving the observed tension in 
the 5 s -meson sector. We emphasize that after minor modifications, our general results can 
be applied to other dimension-six operators involving quarks and leptons. For example, as a 
result of the 90% CL limit B(B + -»■ K + t ± ^) < 7.7 ■ 10" 5 [108], the direct bounds on the 
Wilson coefficients of the set of (sb)(f/j,) operators turn out to be roughly a factor of 7 stronger 
than those in the (sb)(fr) case. Possible effects of (sb)(f/j l ) operators are therefore generically 
too small to lead to a notable improvement of the tension present in the current D0 data. 

Our model-independent study of non-standard effects in r* 2 is complemented by an anal- 
ysis of B s -B s mixing in two explicit standard model extensions. We consider in detail the 
corrections to M( 2 and T\ 2 due to a SU (2) singlet scalar leptoquark and a left-handed Z' boson 
with a mass at or above the electroweak scale. In both cases, we find that depletions of order 
one in AT S and changes by more than two orders of magnitude in a s j s can occur, while a notable 
enhancement of the former quantity with respect to its standard model value is not possible. 
The large effects do, however, not occur simultaneously, since the non-standard contributions 
to M( 2 and Tf 2 are strongly correlated in the considered new-physics scenarios. While Ar s 
becomes maximal (minimal) for 0}/^ = 0° {(j> s ji^ = —180°) and the predictions form a band 
with a cosine-like shape in the ^V^-Ar 8 plane, a s j s is minimal for (fij/^ = 0°, —180° and 
maximal for (fij,^ = —90°, and the accessible parameter space in the ^j/^^^fs plane is 
bounded by sine curves. In turn, neither a SU(2) singlet scalar leptoquark nor a left-handed 
Z' boson renders a good description of the current data. The same conclusion also applies to 
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Operator 


Bound on A (C A = 1) 


Bound on C A (A = 1 TeV) 


Observable 


(sP A b)(fP B r) 


1.3 TeV 


6.6 ■ lO" 1 


B + K+t+t- 


(srP A b)(f^P B r) 


1.0 TeV 


1.1 


B+ K+t+t- 


(sa^P L b)(fa^P L r) 


3.2 TeV 


9.2 • 10~ 2 




(sa^P R b)(fa^P R r) 


2.7 TeV 


1.3 - 10 1 


B^X sl 



Table 3: Model- independent limits on the Wilson coefficients of all possible four- 
fermion operators with flavor structure (s6)(fr). The second column shows the bound 
on the suppression scale A assuming an effective coupling strength of Cf = 1, while the 
third column gives the value for C A for a new-physics scale of A = 1 TeV. Notice that 
only the bounds on the tensor operators depend on the specific chirality A,B = L, R. 



vector-like leptoquarks and Z' bosons with right-handed or both left- and right-handed cou- 
plings. In fact, we have shown that the pattern of deviations found in the case of leptoquarks 
and Z' bosons is a feature of all new-physics model with real (M{ i 2 ) N p/(r^ 2 ) N p and that in 
this class of models the measurement of AM S generically puts stringent constraints on both 
Ar s and a s j s . As it turns out, these bounds are weakest if the ratio r^p = (M[ 2 )np / (X s 12 )np is 
positive and as small as possible. Since on dimensional grounds tnp scales as the square of the 
new-physics scale, this general observation implies that standard model extensions that aim 
at a very good description of the Tevatron data should have new degrees of freedom below 
the electroweak scale and/or be equipped with a (Glashow-Iliopoulos-Maiani-like) mechanism 
that renders the contribution to M\ 2 small. Furthermore, models in which (M^np is gener- 
ated beyond Born level seem more promising, since in such a case tnp is suppressed by a loop 
factor with respect to the case where (M 1 s 2 )np arises already at tree level. 

The above discussion suggests that while the presence of absorptive new-physics allows for 
a better description of the .B^-meson data, a full explanation of the observed discrepancies 
seems not even possible for the most general case (M 12 ) NP 7^ and (r 12 ) NP ^ 0. Numerically, 
we find for example, that the addition of a single (sb)(fr) scalar operator on top of dispersive 
new physics with (Mi 2 )np ^ 0, can only improve the quality of the fit to the latest set of 
measurements (i.e., ([I]), Q, and Q) to x 2 = 1-4 compared to \ 2 = 3-4 within the SM. This 



might indicate that the high central value of A h SL observed by the D0 collaboration is (partly) 
due to a statistical fluctuation. Future improvements in the measurement of the CP phase 
^j/ixj) an d> * n particular, a first determination of the difference a s ^ s — aj s between the B s and 
Bd semileptonic asymmetries by LHCb, are soon expected to shed light on this issue. Will 
these measurements continue to make a case for new physics in rf 2 ? Time will tell! 
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A Operator Mixing 



In the following we present the analytic expressions that are needed to evolve the Wilson 
coefficients from the high-energy scale A down to the low-energy scale m^. In the basis 



(Qs,ALi Qv,ALi Qt,A, Q7,A, Q%,a) 

of these operators, reads 



with A = L, R the LO ADM, which describes the mixing 
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where ft = 11 — 2/3 Nf and Nf denotes the number active quark flavors. Notice that the 
self-mixing of Qs,al is, up to an overall sign, equal to the one of the quark mass. This has 
to be the case, because a composite operator built out of a scalar current and a quark mass 
is conserved and thus has zero anomalous dimension. Furthermore, the self- mixings of Qt,a 
and Q7 t A are the same, up to the factor —2ft (which is due to the normalization of Qt,a)i 



since the tau mass does not run in QCD. Of course, the result (Al), which after a suitable 
replacement of color and charge factors agrees with [531 156"] . also holds for the operators 
{Qs,ARi Qv,ARi Qt,a, Qj,Ai Qq,a) given that gluon interactions conserve chirality. 

Solving the RGEs, we find that the Wilson coefficients at the scale fii < fiQ are given in 
terms of their initial conditions evaluated at /xo by 

C s ,aM=V' / ^C s ,aM, 



Cv,al{^i) — Cv,al{^o) , 
CrM = C 7 , A (»o) + 4- (1 - T 1 ) V m3M C T M , 

Ps 

2 

Cm 0*0 = v' 1 Cm (/■*<)) + ^5- (1 - V 1 ) (C v ,al(Vo) + C V;A rM) , 



(A2) 
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where rj = a s (/xo) / Q! s (a*i)- These results hold if one does not cross a quark threshold in between 
fio and n\. Generalizing them to include threshold effects is however straightforward, since 
the evolution factorizes and finite matching corrections are not present at LO. 



B Kinematic Functions 



In this appendix, we collect the analytic results for the functions Mf s (s) and Mf- + (s), which 
appear in the formulas (30) and (35) for the differential branching ratios of the B — > X s t + t~ 
and B — > K + t + t~ decays. 

In the case of the inclusive b — > st + t~ transition, the relevant kinematic functions are 
given by the following expressions [50] 
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In the case of the exclusive decay, the relevant expressions take instead the form [51] 
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{M B+ + \I K . 



2ml 



1 + — M \(s,M' B+ ,M' K+ )ti(s)g K+ (s) 



with 



g K+ (s) = (3 ^/\(s,M 2 B+ ,M 2 K+ ) , (B5) 

and M K + = 493.677 MeV [35]. The s-dependence of the B + — > K + form factors fi{s) entering 
(B4) is modeled using the results obtained in [109] . Explicitly, we employ (s = s/GeV 2 ) 

f + (s) = (-5.6314 + 6.9718 V33.3258 - s + 0.6206 s) d^s) d 2 (s) , 



2 /If 2 1 



f T (s) = (-7.7322 + 8.2175 V33.3258 - s + 0.7643 s) dx(s) d 2 (s) , (B6) 
f (s) = (-30.040 + 11.200 V33.3258 - s + 1.6269 s) d 2 (s) , 
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withdi(s) = (l-s/Aff.,^)- 1 , d 2 {s) = (4.3173 + V33.3258 - I) 2 , andM B * ( i-) = 5.412 GeV. 



C Matrix Elements for b 



This appendix contains the analytic results for the b — > s£ + £ matrix elements of the complete 
set of (s6)(rr) operators 

Q^isTlMfT^T), (CI) 

with Ti b ® ri r = {P A ® P B , YPa ® o^-Pa ® <V-Pa}- The diagram with the closed tau 

loop shown on the right-hand side in Figure [4] yields the following amplitude for the operator 
insertion 



A = iMQ T (sV i sb b)(I 1 c 



A + InXf, - -M 7 (s,x T ) 



+ 



i (A + lnx M ) - ^M 9 (s,x 7 



Tr [ 7a r^ T ] } Q 



(C2) 



Here C,- is the relevant Wilson coefficient and 



— — ^- t- V t *V t b 



A = - - 7s + In (4tt) 



M 2 ' 



(C3) 



y/2 4vr 

and 7s ~ 0.577216 denotes the Euler-Mascheroni constant. The functions Mj(s,x T ) and 



M g (s,x T ) entering (C2) are given by 

M 7 (s, x T ) = -8 + 4 In x T + 4 
M 7 (0,x r ) = 41nx r , 

and 

M 9 (s,x T ) = 



(C4) 



^(5 + 3y) + hnx T + ^(2 + y)g(y), 



2 10 
M 9 (s,0) = -(lns-m)- — 



(C5) 



where y = 4mJ/s and 



Ky) = v 7 ! 1 -2/1 



In 



2 tan" 1 



1 + ^T^ 
1 



— 27T 



2/ < 1, 
y > 1. 



(C6) 



The terms Tr[7 Q ^r^ T ] and Tr [7 Q r^ r ] in (C2) give rise to the contributions proportional 
to the electromagnetic dipole operators Q 7j a and the vector-like semileptonic operators Qg t A, 
respectively. The relevant non-zero Dirac traces are 



Tr YiaiPA = 2 <?a , Tr [7 a 7 M P A ] = 2g a » , Tr [^(T^Pa] = 4z {q^g av - q v g ai i) • (C7) 
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Finally, the tree- level b — > s£ + £ matrix elements involving the operators (39) read 



(£?7,A)tree 
(Q9,A)tree 



Ttl — 

iN2^{s 1 »iP A b){l 1 ,l) 

iM{s^P A b){l llx i). 



{CI 



By combining (|C2[ ), (C7), and (C8) it is a matter of simple algebra to derive both the expres- 
sions given in (41) as well as the relevant entries of the LO ADM in (Al). 



D Matrix Elements for b — >• S77 

In the following, we present the analytic results for the one-loop matrix elements involving an 
insertion of a (sfe)(rr) operator, which enter the LO prediction for B s — > 77. The relevant 
decay amplitude has the following general structure 



3 



f Bs { A + (2k! ■e 2 k 2 -e 1 - M Bs e t ■ e 2 ) - 2iA_ e{k u k 2 , £1, e 2 ) } . 



(Dl) 



where M has been defined in ( C3 ) , F„ v and F, 



1 2 p x P are the photon field strength 
tensor and its dual, k\ j2 and £^2 denote the momenta and the polarization vectors of the 
two photons, and we have introduced e(a,b,c,d) = e^ v \ p a^b u c x d p . The subscripts ± on the 
coefficients A±, denote the CP properties of the corresponding two-photon final state with A + 
(A_) being proportional to the parallel (perpendicular) spin polarization e\ ■ e 2 (e\ x e 2 ). 

We find that the non-vanishing matrix elements (Qi) = (jj\Qi\B s ) arising from operator 
insertions into the Feynman diagram of Figure El are given by (y T = rn^/M^J 



4y T ) KVt)] (2An ■ £ 2 k 2 - e x - M 2 Bs e x ■ e 2 ) 



=F [2 + h(y T )] ie(ki,k2,e 1 ,e 2 ) 



(D2) 



{Qv,la) 



a 
Arc 



fB s Ql{ ±h(y T )ie(k 1 ,k 2 ,e 1 ,e 2 ) 



with 



h(y) 



1 + y In 2 



(D3) 

L (A — R). The 

expressions (D2) can be shown to resemble the results of [86"| 189]. The formulas that apply in 



Here a 



1/137.036 and the upper (lower) signs hold in the case A 



the case of Qs,ra and Qv,ra are obtained from (D2) by reversing the overall sign. Notice that 
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although the 1PI diagrams with an insertion of a tensor operator with flavor content (s6)(rr) 



evaluate to zero, the operators Qt,a contribute to the amplitude (Dl), because the tensor 
operators lead to corrections AC 7 to the effective Wilson coefficients of Q 7 . We finally 
mention that (Qsab) scales like the mass of the charged lepton in the loop. This implies 
that (Qs,ab) for y — > 0, since h(y) — > —2 in this limit. The matrix elements of scalar 
operators involving electrons or muons are therefore to first approximation equal to zero. In 
the case of tau leptons the chiral suppression is compensated by the large logarithms In y and 
the numerical factors appearing in h(y) = -2+ (-2 In 2 ?/ + 2vr 2 - Ai-nlny) y + 0(y 2 ). 

Combining (Dl ) and (D2) it follows that the new-physics corrections to the coefficients A± 
take the form 



Ms* 
Xb 



(ACf - AC' 7 cS ) 



3 M 



B s 



4 m T 



[2 + (l-4y T )h(y T )]F+ 



^ (ACf + AC'^) ~ I ^ [2 + h(y T )] Fs + | h(y T )E 



(D4) 



v > 



where 



Ci LL ± Cur ~ C, 



J i,RL 



i.RR ■ 



(D5) 



Notice that only the contributions associated with Q 7 contribute to the B s — > 77 amplitudes 
at leading power in AQco/^b, while the effects arising from the (s6)(rr) insertions are relative 
to the leading terms suppressed by one power of \b/Mb s = C^Aqcd/"^)- 

After summing over the photon polarizations the branching ratio of B s — > 77 becomes 



B(B S -> 77) 



A^ 2 Mi Jl r Bs 



144 7T 



( \{A+) m + AA + \ 2 + |(A_)sm + A^-l 2 ) 



(D6) 



where (A±)sm denotes SM contributions to A±. Explicit formulas for the SM coefficients, 
including leading as well as subleading corrections in Aq CD /m b , can be found in [881 189] - 



E Forward- Scattering Amplitudes 

In this appendix, we give some details on the calculation of the contributions from the complete 



set of operators (21) to r* 2 . Since the decay width is related to the absorptive part of the 
forward scattering amplitude, the off-diagonal element of the decay-width matrix may be 
written as 

T S 12 = ^-(B S \T\B S ), (El) 
with the transition operator T given by 

T = 8G 2 F (V*V tb ) 2 J2 Im / d ' x T [CiQi{x) Cj-Qi(O)] J . (E2) 

1 10 



Here T denotes time ordering, the Wilson coefficients Cj are normalized as in (18) and (22). 
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At leading power in AqcdA^& the contribution to T s l2 is found by computing the matrix 
elements of T between quark states. The corresponding Feynman diagram involving an inser- 
tion of a pair of operators (CI) is depicted in Figure [5] Treating the external quarks as 
on-shell and neglecting the strange-quark mass, we obtain from the discontinuity of the graph 



(r; 2 )< 



1 - AX T 

12 



x 



-Mrwri r ] + f^Tr[ri r [r| r ,&]] 



Tr [r frl a n Tla ] 



l + 2x T 
6ml 



(E3) 



with Ar s and j3 T given in and before (55), respectively. Notice that the prefactor f3 T is related 
to the imaginary part of the Passarino-Veltman two-point integral Bo(ml, m 2 rl m 2 ). Explicitly, 
one has 7r/3 T = Im [Bo(ml, m 2 , m^)\. The result (E3) can be shown to agree with the general 
findings of |Iini Ell [T12]. 

The general result (E3) simplifies if one considers the self- interference of operators only. 
Due to the cyclicity of the trace one has Tr [T l fT [r i fT ,$ b )] = for all i = S, AB, V, AB, T, A. 
In the case of i = S, AB, T, A the only non- vanishing Dirac traces read 

Tr [P LiR P LtR ] = 2, Tr [a^P L>R a xp P L , R ] = 2 {g pX g up - g pp g vX =F ie^ xp ) , 



(E4) 



where the Levi-Civita Tensor e pvXp is related to 75 via 75 = — i/4! e pvXp j^ulxlp- For i = V, AB 
one needs instead 



Tr[ 7 p P L , Rl a YPL,R7 a } = -4g 



[IV 



Tr [YPl,rAYPl,rA] = 2 (2|« ~ m\gT) ■ (E5) 



In order to derive the expressions given in (54), one still has to take into account that 

ie pvXp (sa^P L>R b) (sa Xp P LtR b) = T2 (sa^P L , R b) (sa, u P L>R b) , 
and apply the Fierz identity 



(E6) 



[sa^P L>R b) (sa^P^R b) = {sP L)R b) (sP L>R b) + 2 (s a P L>R bp) (s fi P L , R b 



(E7) 



F Bounds on Ar s and a s ^ s from Ail/f 



In this appendix, we show that in the class of SM extensions with real (M^ 2 ) NP /(r^ 2 ) N p, a 
bound on AM S necessarily results in a limit on both Ar s and a s j s . We start our derivation 
by noting that under the assumption that arg ((Aff 2 )Np/(r^ 2 )Np) is either equal to 0° or 180°, 



the parameters Rm,t introduced in (14) can be expressed in terms of the three ratios 



(M; 



12/NP 



(M- 



r SM 



(M; 



12JSM 



12JSM 



(r 



?~NP 



(AT 



12JNP 



12JSM 



(r 



12/NP 



(Fl) 



and 0g M , as follows 



R 



Al 



\l + r M \ 



l-sgn(r NP )P^r M e^M 
r N p 



1 - sgn (r NP ) | f r M 

r N p 



(F2) 
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Because 0| M is only a fraction of 1°, the approximation made to obtain the final expression 
for i?r is sufficient for all practical purposes. 

In order to simplify the further discussion, we consider the two possible signs of tnp sepa- 
rately. We begin with sgn (tnp) = +!■ From (F2) it is readily seen that in this case, the value 



of Rr becomes maximal (minimal) if tm is real, negative (positive), and as large as possible 
in magnitude. A two-sided bound on AM S or equivalently Rm, 



r>mm ~ z~> , pmax 



(F3) 



with R M m < 1 and -R^ ax > 1 18 thus constrains the possible values of Rr to lie in the interval 
max ( 0, 1 - (R^r ~ l) < Rr < I + (#m x + l) 

V l r Np|y 



\rsM\ 



(F4) 



According to (15) a bound on Rr also restricts the allowed range for Ar s . Since one has 
4>m = 180° and <pr = 0° if Rr becomes maximal, the upper limit in (F4) leads to a lower 
bound on Ar s . On the other hand, the upper limit on Ar s arises if <$m 
tm = Rti n - 1 < -R™ ax - 1- One therefore has 



0° and 



M 



1 + {R 



max 
M 



1) 



|tsm| 
IrNpl 



(Ar s ) SM < (Ar s ) NP < 



1 - (R™ ~ 1) 



Fsm| 
k N p| 



(Ar s ) SM • (F5) 

Given that both Rm and Rr are bounded it is not surprising that also a^ s satisfies a dou- 
ble inequality. Because the expression (15) for a s j s is inversely proportional to Rm, the CP 
asymmetry becomes extremal for Rm = R M m - It is then easy to convince oneself, that for 
this value of Rm the combination of Rr sin ((pM ~ <Pr) is maximized/minimized if tm is purely 
imaginary, corresponding to 4>m — ±90°. After some algebraic simplifications, one arrives at 



pmin 

M 



1 + 



Fsm| 



la's) 



fsJSM 



pmin 

M 



1 + 



k N p| 



l /JSM 
^SM 



(F6) 



Using the same line of reasoning as above, it is also not too difficult to derive the constraints 
for sgn (rNp) = —1 that follow from (F3) with R M m < 1 and R™ x > 1. One has to differentiate 



between the case |rsM|/|?"Np| < 1 an d |^sm|/|^np| > 1- in the former case, we obtain 
max ( 0, 1 - (Rfr + l) P4] < Rr < I + (HST - l) ^ ' 



-1 + (R™ + 1) 



|tsm| 
k N p| 



(Ar s ) SM < (Ar s ) NP < 



i + (Rtr - 1) 



Tnp 

I^smI 
k N p| 



(F7) 



(Ar.) 



SM 



while in the latter case the relations 



max ( 0, 1 + (R™ D - l) j^J ] < R r < -1 + (i^ ax + l) 

V l r Np|/ 



|tsm| 
Insrpl 



18 This is the relevant case in view of the good agreement between (AM s )sm and measured value of the mass 
difference. It is straightforward to extend the given formulas to the other possible cases. 
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1 + (R™ n - 1) f^i 



(Ar s ) SM < (Ar s ) NP < 



•i + (i?- ax + i) 



\ r sm| 
|r N p| 



(AT.) 



SM ■ 



(F8) 



apply. The constraint on the CP asymmetry is independent of whether |7"sm|/|?"np| is smaller 
or bigger than 1. It takes the simple form 



pmm 



\rsu\ 



r N p 



(fl/JSM , s . 1 

— Q < KJnp < 



y SM 



pmin 



r N p 



i/JSM 



(F9) 



The formulas (F4) to (F9) have a couple of features that are worth mentioning. First, the 



possible variations increase with decreasing IrsMl/I^Npl- Second, the limits on Ar s and a^ s 
are less stringent if sgn (r^p) = +1. The only exception is the upper bound on Ar s , which 
becomes weakest for sgn(r NP ) = —1 and |r SM |/|rNp| > 1. In new-physics scenarios with real 
(Mf 2 ) NP /(r^ 2 ) NP , non-standard effects in AT S and a s j s are hence the least constrained by the 
measurement of AM S , if the ratio tsm/^np is positive and as small as possible. 
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